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Abstract 

Theories  of  Alfven  resonance  and  Kelvin-Helmholtz  instability  are  reviewed  to 
elucidate  a  possible  mechanism  of  the  long-period  continuous  magnetic  pulsations 
(pc3  ~  pc5). 
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1.      Introduction 

It  is  now  known  that  the  earth's  magnetic  field  fluctuates  with  a  wide  range  of  frequencies. 
Aside  from  extremely  long  time  scales  associated  with  the  earth's  interior  physics  (such  as 
dipole  reversal  of  the  earth's  magnetic  field),  the  time  scale  of  the  magnetic  field  fluctuations 
ranges  from  several  days,  caused  by  magnetic  storms,  to  less  than  a  second,  caused  by  Larmor 
motion  of  protons.  The  temporal  changes  in  the  geomagnetic  field  with  the  range  of  periods 
from  fractions  of  a  second  to  several  minutes,  with  which  we  are  particularly  concerned 
in  this  review,  are  called  magnetic  pulsations.  These  pulsations  have  been  measured  with 
magnetometers  placed  on  the  earth's  surface  (ground  observation)  and  on  space-craft  (in-situ 
observation). 

The  magnetic  pulsations  have  been  known  for  over  a  hundred  years.  Historically  such 
rapid  changes  in  the  earth's  magnetic  field  were  first  reported  by  B.  Stewart  in  1861.  He 
reported  a  large  magnetic  storm  recorded  at  the  Kew  Observatory  in  England.  It  is  shown 
in  this  report  that  the  magnitude  of  the  magnetic  field  fluctuations  is  of  several  hundred  7 
(I7  =  10~^gauss)  and  their  typical  time  scale  is  on  the  order  of  few  minutes.  In  particular, 
during  the  International  Geophysical  Year  (IGY)  in  1957-1958  and  in  the  succeeding  years, 
the  international  collaboration  on  the  earth's  geophysical  study  became  established  and  the 
global  observation  of  the  geomagnetic  field  variations  such  as  magnetic  pulsations  began  to 
be  made. 

Based  on  such  global  observation,  the  earth's  magnetic  field  variations  have  been  classi- 
fied into  categories  according  to  the  frequency.  One  of  the  simplest  and  most  widely  used 
classifications  is  based  on  nine  categories  according  to  the  regularity  and  the  period.  In 
this  classification  the  pulsations  pel  ~  pc6,  which  have  rather  well-defined  spectral  peaks, 
are  called  continuous  pulsations  while  the  pulsations  pil  ~  pi3,  which  have  a  wide  spectral 


range,  are  called  irregular  pulsations.  This  is  summarized  in  Table  1,  where  pc6  and  pi3 
are  not  referred  to  but  another  kind  of  classification  sc  and  si  is  used.  A  more  detailed 
classification  has  been  proposed  by  Saito  [1976],  where  time  evolution  of  frequency  of  each 
pulsation  is  also  taken  into  account  (Table  2).  In  the  short  period  range  for  pel  and  pil,  the 
subclassification  based  on  the  time  evolution  of  frequencies  is  considered  to  be  particularly 
important. 

The  detailed  mechanism  of  the  variations  of  the  earth's  magnetic  field  has  not  been  under- 
stood completely  as  yet.  However,  it  is  widely  beheved  that  the  pc  pulsations  (in  particular 
pc3  ~  pc5)  are  caused  by  the  resonant  oscillations  of  the  earth's  magnetic  field  lines  since 
they  have  distinctive  periodicity.  Because  of  the  high  conductivity  of  the  ionosphere,  the 
shear  ALfven  wave  travehng  along  the  magnetic  field  hne  is  expected  to  be  reflected  at  the 
intersection  of  the  field  line  and  the  ionosphere  and  a  standing  wave  will  be  excited  on  the 
field  line  (field  line  resonance).  Therefore,  the  fundamental  period  of  the  oscillation  is  given 
by 

r  =  2/^,  (1.1) 

where  v^  denotes  the  Alfven  velocity  and  the  integration  is  to  be  taken  between  the  two 
conjugate  points  along  a  given  field  fine.  Therefore  the  period  given  by  Eq.  (1.1)  depends 
on  geomagnetic  latitude  or  the  L-value  of  the  field  line.  Figure  1  shows  how  the  period  is 
expected  to  vary  with  latitude  if  Eq.  (1.1)  holds  for  every  field  line  [Saito,  1976]. 

This  simple  theory  of  the  field  line  resonance,  however,  does  not  bear  out  the  observation 
that  the  pulsation  with  the  same  frequency  is  often  recorded  simultaneously  over  a  latitude 
range  of  several  degrees  [Hirasawa,  1970;  Samson  and  Rostoker,  1972].  Indeed,  as  we  will 
show  in  Sec.  2,  field  fines  having  different  eigenfrequencies  can  be  coupled  together  and 
oscillate  as  a  mass.  Therefore  the  period  given  by  Eq.  (1.1)  does  not  necessarily  indicate 
the  observed  period  of  magnetic  pulsations  and  some  other  mechanism  may  be  responsible 
to  make  certain  spectral  peaks  stand  out  in  the  pc  spectrum.  It  is  now  believed  that  there 


are  basically  two  types  of  mechanisms  causing  spectral  peaks  in  the  pc  oscillations.  One 
mechanism  is  an  external  one,  namely  a  monochromatic  oscillation  is  excited  away  from  the 
resonant  field  Hne,  probably  outside  or  at  the  boundary  of  the  magnetosphere,  and  then  it  is 
transmitted  toward  the  earth  as  a  compressional  wave  and  resonates  on  a  field  line.  One  of  the 
most  prominent  candidates  for  the  external  wave  source  is  the  Kelvin- Helmholtz  instability, 
which  is  expected  to  occur  at  the  magnetopause  or  the  boundary  between  the  solar  wind  in 
the  magnetosheath  and  the  plasma  in  the  magnetosphere.  Once  such  a  monochromatic  wave 
propagates  in  the  magnetosphere,  Alfven  waves  are  excited  in  response  to  the  source  wave 
in  the  region  where  the  eigenperiod  of  the  geomagnetic  field  line  coincides  with  the  period  of 
the  source  wave.  This  mechanism  is  called  Alfven  resonance  and  due  to  the  nonuniformity 
of  the  magnetic  fields.  This  phenomenon  heis  been  applied  to  many  situations  in  laboratory 
and  space  and  astrophysical  plaismas  over  the  past  twenty  years  [Tataronis  and  Grossmann, 
1973;  Grossmann  and  Tataronis,  1973;  and  Hameiri  and  Hammer,  1979]. 

The  other  mechanism  which  might  be  responsible  for  the  spectral  peaks  in  the  magnetic 
pulsations  is  an  internal  one,  namely  excitation  of  a  monochromatic  wave  occurs  inside 
the  magnetosphere.  One  of  the  widely  accepted  theories  of  this  mechanism  is  excitation 
of  the  surface  wave  on  the  surface  where  the  Alfven  velocity  changes  abruptly  [Chen  and 
HcLsegawa,  1974b].  This  mechanism  explains  plasma-pause  associated  magnetic  pulsations 
(the  frequency  of  which  falls  into  pc3  ~  pc5),  which  are  weakly  damped  coherent  oscillations 
at  L  ~  3  and  are  accompanied  by  impulse-like  perturbations  at  Z  ~  4  [Lanzerotti  et  al., 
1973]. 

In  this  review,  we  focus  on  the  first  mechanism  of  the  pc-pulsations  mentioned  above 
[Nishida,  1978;  Lanzerotti  and  Southwood,  1979;  Rostoker,  1979;  Haerendel  and  Paschmann, 
1982;  Southwood  and  Hughes,  1983];  in  particular,  the  Kelvin-Helmholtz  instability  is  con- 
sidered as  a  possible  wave  source.  The  rest  of  this  paper  is  organized  as  follows.  In  Sec.  2 
the  theory  of  Alfven  resonance  is  presented.    In  Sec.  3  the  magnetohydrodynamic  (MHD) 


Kelvin-Helmholtz  instability  is  discussed.  The  last  section  contains  summary  and  discus- 
sion. 

2.      Resonant  Excitation  of  Alfven  Waves 

We  now  consider  the  response  of  the  magnetosphere  to  a  monochromatic  source  wave,  assum- 
ing that  such  a  source  wave  may  be  generated  by  an  instability  such  as  the  Kelvin-Helmholtz 
instability  at  the  magnetopause.  A  coupling  between  an  active  process  releasing  free  energy 
such  as  the  Kelvin-Helmholtz  instability  and  a  passive  process  such  as  the  resonant  excita- 
tion of  an  Alfven  wave  (which  is  observed  as  geomagnetic  pulsations)  was  first  studied  by 
Southwood  [1974],  using  a  straight  field  line  model,  and  also  independently  by  Chen  and 
Hasegawa  [1974a],  using  the  dipole  field  model.  In  this  section,  we  mainly  follow  the  work 
of  Chen  and  Hasegawa  [1974a]  (see  also  Hasegawa  and  Chen  [1974]). 

We  assume  that  the  magnetosphere  can  be  described  by  the  ideal  MHD  equations: 

d 

p—v  =  -Vp  -I-  ( V  X  B)  X  B,  (2.1a) 

at 


(9R 

—  =  V  X  (V  X  B),  (2.1b) 


where 


and 


^  +  7pV  .  V  =  0,  (2.1c) 


^  =  |  +  v.V,  (2.1d) 


VB  =  0,  (2.1e) 


p  =  Ap\  (2.1f) 


Here  p,  p,  B,  and  v  denote  the  mass  density,  the  pressure,  the  magnetic  field,  and  the  velocity- 
field,  respectively,  7  denotes  the  ratio  of  the  specific  heats,  and  A  is  a  constant.  We  write  each 
dependent  variable  as  a  sum  of  the  time-independent  mean  part,  denoted  by  subscript  0,  and 
the  small  fluctuating  part,  denoted  by  subscript  1.  For  example,  B  =  Bq  +  Bi.  Assuming 
that  the  mean  fluid  velocity  vanishes,  we  write  v(=  Vi)  as  the  fluctuating  fluid  velocity. 
Since  the  following  hnearized  equations  do  not  involve  the  fluctuating  mass  density,  we  write 
P  (=  Po)  as  the  mean  part  of  the  mass  density. 

Assuming  that  the  mean  quantities  Bo,  po,  and  p  are  zeroth  order  quantities  while  the 
fluctuating  quantities  are  first  order  quantities,  we  linearize  the  system  (2.1): 

pi    =     -V(pi  +  Bo-Bi)  +  (Bi-V)Bo  +  (Bo-V)Bi  (2.2a) 

Bi     =     -Bo(V-.0  +  (Bo-V)C-(^-V)Bo  (2.2b) 

p:     =     -(^.V)po-7Po(V-0,  (2.2c) 

where  ^  denotes  the  plasma  displacement,  defined  by  ^  =  v  and  the  equilibrium  condition 
Vpo  =  (V  X  Bo)  X  Bq  is  used.  Substituting  Eq.  (2.2b)  into  Eq.  (2.2a),  we  obtain 

p'k  -  (Bo  •  V)'  e  =  -V  (pi  +  Bo  ■  Bi )  +  C,  (2.3) 

where 

C  =  -  (Bo  •  V)  {{i  ■  V)Bo  +  Bo(  V  •  ^)}  +  (Bi  •  V)  Bo-  (2.4) 

We  note  here  that  the  dispersion  relation  of  the  shear  Alfven  wave  is  given  by 

p$^-(Bo-V)^e^     =    0 

and  the  dispersion  relation  of  the  surface  wave  excited  by  the  Kelvin-Helmholtz  instability 
in  an  incompressible  pleisma  is  given  by  [Sen,  1963] 

V2(pi-fBo-Bi)    =    0, 


which  will  be  shown  in  Sec.  3  [Eq.  (3.5)].  Thus  Eq.  (2.3)  represents  a  coupHng  between 
the  shear  Alfven  wave  and  the  surface  Alfven  wave  with  the  strength  of  the  coupling  being 
determined  by  C. 

In  order  to  solve  the  system  of  equations  (2.2c),  (2.3),  and  (2.4)  in  a  special  case,  we  now 
assume  a  simple  geometry.  Suppose  all  mean  quantities  depend  only  upon  the  y-coordinate, 
the  straight  magnetic  field  hne  is  in  the  z-direction  and  all  fluctuating  quantities  depend  on 
time  and  space  in  a  form  f{y)  exp  i  (k^x  +  k\\z  —  ujij ,  we  obtain 

from  Eqs.  (2.2c),  (2.3),  and  (2.4).  In  this  geometry,  the  y-direction  represents  the  radially 
outward  direction  in  the  actual  magnetosphere,  the  2-direction  represents  the  south-north 
direction  and  the  x-direction  represents  the  east-west  direction.  In  Eq.  (2.5),  ^y  is  the 
j/-component  of  ^, 

e{y)    =    u'p{y)-klB'{y), 

a{y)     =     1  +  7/3  + 


u^p  -  fl3k^.B^ 


2D2' 


anc 


/3{y)  =   ^, 


B    =    |Bo|. 

Furthermore,  we  assume  a  wave  locahzed  in  the  east-west  direction,  i.e.,  k±  >  A:||  ~  u/va, 
where  v^  =  B/ y/p  is  the  Alfven  velocity.  Then  Eq.  (2.5)  becomes 

^  +  ^  ^  _  kl^y  =  0,  (2.6) 

dy^         dy     dy 

where  we  used  aB'^k\  >  |£:|.    For  a  given  source  wave  frequMicy  u  (which  generally  has  a 
small  positive  imaginary  part  indicating  damping)  the  position  t/o  of  the  resonant  field  line 
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is  given  by  the  condition  Ree  (t/o)  =  0.  Away  from  the  resonant  field  line,  the  plasma  and 
the  magnetic  field  are  almost  homogeneous  (i.e.,  dlne/dy  ~  0)  except  for  some  discontinuity 
surface  such  as  the  magnetopause.  In  such  region,  Eq.  (2.6)  is  of  the  form  V^i^j,  =  d'^^y/dy^  — 
^"liy  —  Oi  which  gives  the  surface  wave  solution  ^j,  =  ^y  exp  (—  \k±_\  [y  —  yd))-  Here  yd  denotes 
the  position  of  the  discontinuity  surface.  The  strong  coupling  of  the  surface  wave  and  the 
local  Alfven  wave  at  y  =  yo  occurs  when  the  discontinuity  surface  is  close  to  the  resonant 
surface. 

Near  the  resonant  field  fine,  we  may  obtain  the  solution  of  Eq.  (2.6)  by  expanding  e{y) 
around  the  resonant  point  yo 

de 


e(y)  =  £(yo)  + 


dy 


{y  -  yo) 


y=vo 


If  weak  dissipation  is  present,  then  the  source  wave  has  a  small  positive  imaginary  part 
Imu;(<  Reu).  In  this  case,  \{s:/de/dy)y^y^\  ~  |(£t75£r/52/)y=„J  is  small  compared  with  a 
typical  scale  length  of  the  variation  of  plasma  parameters.  Here  e^  =  Ree  and  e,  =  Ime. 
Using  the  following  parameter 


7'  = 


£, 


der/dy 


y=jra 


we  therefore  rewrite  Eq.  (2.5)  near  yo  as 


<P^y 


+ 


_^_fc2^     -0 


(2.7) 


dy^       y  -  2/0  +  J7'  dy 

The  Alfven  speed  va  =  B/y/p  is  believed  to  be  an  increasing  function  of  the  distance  toward 
the  earth  (i.e.,  y  in  our  case)  except  at  the  plasma  pause.  Therefore,  the  y-derivative  of  e^ 
at  y  =  yo 


dCr 
dy 


y=yo 


_        \  dp  d_ 

~      ,P  dy  ""        dy 


—  Km 


y=yo 


=   -4("J^-n) 


dy 


y=yo 


is  negative  except  at  the  plasma  pause,  and  so  is  the  parameter  7'.  We  note  that  k\\  is 
generally  an  increasing  function  of  the  distance  toward  the  earth,  although  in  our  planar 
geometry  k\\  is  a  constant. 

The  solution  of  Eq.  (2.7)  is  given  by 

(y  =  do  {k^  {y-yo  +  h'))  +  DKo  {k^  {y-yo  +  n)) , 

where  /q  a^nd  A'o  are  modified  Bessel  functions  and  C  and  D  are  constants.  The  corresponding 
eigenmode  is  shown  to  be  continuous.  With  a  monochromatic  wave  source,  the  solution  of 
Eq.  (2.7)  near  y  =  yo  is  approximately  given  by 

^y    ~    Cln(y-l  +  i7?),  (2.8) 

where  Y  =  y/yo  and  77  =  -y'/yo-  It  can  be  shown  that  the  wave  near  the  resonant  field  line 
y  =  yo  is  noncompressional  [Hasegawa  and  Chen,  1974],  a^  one  would  expect  for  a  shear 
ALfven  mode.  Therefore  the  equation  V  •  ^  =  0,  together  with  the  assumption  k±  ^  k\\, 
yields  the  expression  of  the  x-component  of  (  near  j/o 

^'    -     klioY-l  +  ir,-  ^^-^^ 

We  now  discuss  the  polarization  of  the  wave  obtained  above.  The  relationship  between 
^  and  Bi  is  given  from  Eq.  (2.2b) 

{Bix  =  ik\\Bo^x 
Biy  =  ik\\Boiy, 

where  the  subscripts  x  and  y  indicate  the  x  and  y  components  of  the  associated  quantities, 

respectively.  Therefore,  near  the  resonant  field  hne,  we  obtain  from  Eqs.  (2.8)  and  (2.9) 

rf\y 

=    g  +  ih,  (2.10) 

where 


9    =     k^yo  [l\n  {ay' +  rj')+ AY t^n-'^y 
h     =    k^y,(rjte.n-'^-^ln{AY'  +  rj')y 


and 

AY    =    Y-1. 

Since  Ay  ~  0  near  the  resonant  field  line,  we  have 

9     -     k_iyo-\D.Ti^ 
h     ~    k_Lyo\T]\- sign  AY. 
The  signs  of  the  functions  g  and  h  are  then  given  by 

signi?    =     -sign(tx'?)  (2.11) 

sign/i    =    s\gn{k^AY)  (2.12) 

since  r/^  <C  1  or  hi;?^  <  0.  From  Eq.  (2.10)  we  have 

p     J"^  =  5  -  /i  tan  [kxx  +  k\\z  -  {Reu)tj  , 

so  the  sense  of  polarization  is  determined  by  the  sign  of  h.  From  Eq.  (2.12),  it  is  clear  that 
the  sense  of  polarization  switches  across  the  resonant  field  line  at  y  =  yo- 

The  ellipticity  and  the  angle  of  polarization  can  be  seen  from  the  equation 

{ReB^yf  -  2g  (Re^iJ  (ReBiJ  +  (5'  +  h'')  (Re^i^)'  =  |Bi^f  h\  (2.13) 

which  is  also  derived  from  Eq.  (2.10).  It  is  easy  to  show  that  the  orbit  of  the  point  (ReBii, 
Refiiy)  given  by  Eq.  (2.13)  is  an  ellipse.  Solving  Eq.  (2.13)  for  ReBi^,  we  obtain 

Re  6,  =  g{KeB,,)  ±  hyj  \B,,\'  -  (Re^i,)'.  (2.14) 

10 


Hence,  as  one  approaches  the  resonant  field  hne  (i.e.  h  —^  0),  the  elhpticity  becomes  large 
and  the  polarization  becomes  linear  at  the  resonant  field  line  {h  =  0).  Figure  2  summarizes 
the  angle  of  polarization  obtained  from  Eq.  (2.14).  In  this  figure,  k±  is  written  as  m,  x  and 
y  are  written  as  4>  and  i),  respectively.  Assuming  that  the  monochromatic  source  wave  is 
generated  by  the  Kelvin-Helmholtz  instabihty  at  the  magnetopause,  the  sign  of  m  changes 
across  the  local  noon:  m  is  negative  (positive)  in  the  morning  (afternoon)  side.  As  shown  in 
Fig.  2,  the  angle  of  the  polarization  switches  across  the  local  noon  since  the  sign  g  changes 
as  the  sign  of  m  changes.  The  amplitude,  the  elhpticity,  and  the  sense  of  the  polarization 
of  the  wave  on  the  morning  side  (m  <  0)  are  shown  in  Fig.  3.  In  this  figure,  the  parameter 
i/{=  y)  corresponds  to  Z,- value,  so  the  horizontal  axis  also  can  be  viewed  as  latitude.  The 
change  of  the  sense  of  the  polarization  given  in  Fig.  3  bears  out  the  observation  shown  in 
Fig.  4  [Samson  et  al.,  1971]. 

3.      Kelvin-Helmholtz  Instability 

In  this  section  we  discuss  the  Kelvin-Helmholtz  instability  at  the  magnetopause,  which  is 
considered  to  be  a  source  of  the  monochromatic  wave  that  leads  to  resonant  oscillations  of 
the  geomagnetic  field  lines.  For  simphcity  we  first  ignore  the  width  of  the  magnetopause. 
It  is  also  assumed  that  there  is  a  uniform  infinite  flow  of  the  magnetized  plasma  on  the 
magnetosheath  side  (Side  I),  whose  velocity  Vq  is  tangential  to  the  magnetopause,  and  the 
plasma  on  the  magnetosphere  side  (Side  II)  is  at  rest.  Giving  a  small  perturbation  to  this 
magnetosheath-magnetosphere  boundary  model,  we  will  follow  the  time  evolution  of  the 
perturbation  [Sen,  1963]. 

We  use  the  incompressible  ideal  MHD  equations,  i.e.,  Eqs.  (2.1a),  (2.1b),  (2. Id),  (2.1e), 
and 

V-v  =  0,  (3.1) 
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together  with  the  boundary  conditions 


V„-Vs      =      0 


[5„]    =    0 
=    0. 


(3.2a) 
(3.2b) 
(3.2c) 


Here  the  subscripts  n  and  t  denote  the  directions  normal  and  tangential  to  the  discontinuity 
surface  or  the  magnetopause,  respectively,  Vg  refers  to  the  (normal)  velocity  of  the  discon- 
tinuity surface  and  the  square  bracket  is  the  change  in  the  value  of  the  enclosed  expression 
across  the  discontinuity.  It  should  be  noted  that  we  obtain  the  boundary  conditions  (3.2a), 
(3.2b),  and  (3.2c)  by  requiring  that  Eq.  (3.1),  V  •  B  =  0  and  the  normal  component  of  the 
momentum  Eq.  (2.1a)  hold  across  the  boundary.  As  in  Sec.  2,  we  write  each  dependent 
variable  as  a  sum  of  the  time-independent  mean  part  and  the  small  fluctuating  part  such  as 
B  =  Bo  +  Bi,  V  =  Vq  +  Vj,  and  p  =  Po  +  Pi,  and  then  linearize  the  system  of  the  equations. 
Here  we  assume  that  po  =  0  on  both  sides  and  Vq  =  0  in  the  magnetosphere  (Side  II). 
In  order  to  distinguish  the  variables  on  Side  I  and  II,  we  use  superscripts  (I)  and  (II);  for 
example,  Bq  denotes  the  mean  magnetic  field  on  Side  I  and  p\  denotes  the  perturbed 
pressure  on  Side  II.  We  choose  a  coordinate  system  in  such  a  way  that  the  y-direction  is 
normal  to  the  discontinuity  surface  (from  the  magnetosphere  (Side  II)  to  the  magnetosheath 
(Side  I)),  the  i-direction  is  in  the  direction  of  the  mean  flow  vj,  in  the  magnetosheath  and 
the  2-direction  is  perpendicular  to  these  two  directions  (Fig.  5).  By  linearizing  the  system 
of  the  equations,  we  obtain 

^'•^  (I  +  v("  .  v)  vl')    =     -Vpl')  +  (V  X  B';-^)  X  B<'\  (3.3a) 


^    =     -(v?.v)b?  +  (b?.v)vS'\ 


dt 
V-vP    =    0, 


(3.3b) 
(3.3c) 
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.(••) 


V-Bl''    =    0,  (3.3d) 

for  both  fluids  i  =  I  and  II.  As  we  have  assumed,  Vq  and  Bq  are  constant  vectors  and 
vl,"^  =  0. 

Rewriting  the  last  term  of  the  right-hand  side  of  Eq.  (3.3a)  as 

(V  X  B<'')  X  B«  =  -V  (bI,''  •  Bi*^)  +  (Bo  •  V)  B^  (3.4) 

and  taking  the  divergence  of  both  sides  of  Eq.  (3.3a),  we  obtain 

V2(p(')  +  Bi'^-BP)=0.  (3.5) 

Applying  Bq  to  both  Eq.  (3.3a)  and  Eq.  (3.3b)  and  ehminating  Bo  •  Vi  from  the  resulting 
equations  yields 

'''  (I  -^ "''  •  ^)'  (^°'  •  ^'') + (^''  •  ^y^'' = °-       ^'-'^ 

We  now  denote  the  x-,  j/-,  and  z-components  of  Bj'^  by  Bjj,  B[y,  and  B^^ ,  respectively. 
Taking  the  y-component  of  Eq.  (3.3a)  and  applying  Bq^  •  V  to  the  resulting  equation,  we 
obtain 

/'(|  +  v^".v)(B(".v).l;>    =    -^  (B^"  .  V)  (p(')  +  B^')  .  Bl") 

+  (B«*)-V)'Bg.  (3.7) 

Eliminating  (bJ,'^  •  V)  vJ'J  from  the  y-component  of  Eq.  (3.3b)  and  Eq.  (3.7)  yields 

Integrating  Eq.  (3.3b)  for  i  =  II  in  time  at  y  =  0,  we  obtain 

Bi"^  -  (bS°^  •  V)  ^    =    0,  (3.9) 

where  ^  denotes  the  infinitesimal  displacement  of  the  discontinuity  surface  (Fig.  5a).  From 
the  boundary  conditions  (3.2a)  and  (3.2b),  the  similar  equation  holds  for  Side  I  at  y  =  0; 
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B?  -  (b?  •  v)  ^    =    0.  (3.10) 

From  Eq.  (3.5),  we  also  have 

[p:  +  Bo  •  Bi]  =  0.  (3.11) 

We  now  consider  the  normal  mode  solutions  of  this  linearized  system,  assuming  that  all 
dependent  variables  have  the  following  dependence  on  i,  r  =  {x,z)  and  y: 

exp  i  {ut  —  (k  •  r)  —  7ij/)    in  fluid  (I)      {y  >  0) 

expi{ujt-{k-r)  +  -f2y)    in  fluid  (II)    (y  <  0),  (3.12) 

where  k  is  the  propagation  vector,  u;  is  the  complex  frequency,  and  71  and  72  are  real  positive 
constants  so  that  the  solutions  decay  as  |j/|  — >  00.  The  propagation  vector  k  is  taken  to  be 
the  same  in  both  fluids.  From  Eq.  (3.5),  we  obtain 

71  =  72  =  t  =  ik|.  (3.13) 

Using  the  dependence  given  by  Eqs.  (3.12)  and  (3.13),  we  derive  the  following  Unear  algebraic 
homogeneous  equations  in  the  matrix  form: 
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0 
:(k.B<')) 


-ik 


(k.B<"))^ 

0 

(k.BS,"') 


0 
0 

-1 


-p(')(u.-k-vo)^ 
0 

.>(kB(')) 

0 
0 
0 

1 


-p(')(^-k.Vo)^+(k.B^"y 


-,fc(k.B^"') 


0 
0 

-1 


-p(")u,2+(k.B<"')' 


tkB 


(I) 


.kB 


(II) 


Bj.".B(/' 
3(11). B(ii) 

r(I) 

h(ii) 


=  0 


(3.14) 


Here  Vq  =  Vq    and  Vq     —  0.  The  condition  that  nontrivial  solutions  exist  yields  the  following 
dispersion  relations: 


u    =     ±A:Cl")cos(<^„-V), 


u 


Vj  cos  xl^±Ci^  cos  {<j)l-Tp)] 


and 


u 


=     k- 


Pi 


±    i- 


Pi  +  P2 
k 


vq  cos  xp 


(3.15a) 
(3.15b) 

(3.15c) 


Pi  +  P2 


{p,P2vl  COS^  0  -  (/9i  +  P2)  {p,C^^^'  cos^  (<^i  -  0)  +  ^2Ci°)^  cos^  (^n  -  V-))} 


1/2 


where  C^'^  =  |b[,''|  /^pW  (j  =  I,II)  is  the  Alfven  speed  of  the  fluid  i,  uq  =  |vo^^|,  (f>,  is  the 
angle  between  x-axis  and  the  magnetic  field  B^'^  [i  =  I,  II)  and  ip  is  the  angle  between  x-axis 
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and  the  propagation  vector  k  (Fig.  5).  Clearly  the  dispersion  relations  (3.15a)  and  (3.15b) 
represent  stable  solutions.  From  the  dispersion  relation  (3.15c),  we  obtain  the  following 
condition  of  instabihty: 


Pi  +  P2 
which  can  also  be  written  as 


vlcos^^;    >     Pl±Pl[B^Jhcos^,f>l-^P)  +  Bl|'^'cos^<Pu-^P)),  (3.16a) 

P1P2  ' 


<"-»)'    >    ^75^(''""('<-Ci")'  +  /'™(k-C™)'),  (3.16b) 

where  C^,''  =  "B^q^  j -J^)  {i  =  I, II).  When  this  condition  is  met,  the  instabihty  sets  in  at 
the  magnetopause.  This  condition  tends  to  be  satisfied  near  the  equatorial  plane  of  the 
dayside,  where  k  is  parallel  to  the  streaming  velocity  Vq  and  perpendicular  to  the  magnetic 
field  B^"'  in  the  magnetosphere,  i.e.,  0  =  0,  (^i  =  7r/2.  By  taking  into  account  the  fact  that 
p^  '  ^  p^  \  the  foregoing  condition  for  the  Kelvin-Helmholtz  instabihty  with  ip  =  0  and 
'i>u  =  ""/S  becomes 

vl>B^^^'cos^(i>i/p^^^  (3.17) 

and  the  growth  rate  Imo;  of  the  instability  is  given  by 

We  note  here  that  the  right-hand  side  of  Eq.  (3.17)  represents  the  stabilizing  effect  of  the 
magnetosheath  magnetic  field. 

In  the  previous  model,  we  did  not  take  into  account  effects  of  finite  thickness  of  the 
magnetopause.  However,  as  shown  below,  waves  with  wavelengths  comparable  to  the  thick- 
ness of  the  magnetopause  are  significantly  affected  by  finite  thickness  of  the  magnetopause. 
For  example,  the  hnear  growth  rate  shown  in  Eq.  (3.18)  shows  that  a  perturbation  with 
shorter  wavelength  (larger  k)  has  a  larger  growth  rate.  This  apparent  counter-intuitive  large 
growth  rate  of  a  wave  with  large  k  is  reduced  significantly  by  the  sheared  flow  inside  the 
magnetopause. 
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My)  = 


Ong  and  Roderick  [1972]  used  the  following  model  of  a  shear  layer  of  thickness  2c?  in  the 
x-z  plane  as  shown  in  Fig.  6.  This  configuration  models  the  flow  field  and  magnetic  field  on 
the  dusk  side  in  the  equatorial  plane;  regime  (I)  and  (III)  represent  the  magnetosheath  and 
the  magnetopause,  respectively,  with  constant  magnetic  fields  Bq  and  Bq  and  constant 
flow  field  ±V/2.  The  equihbrium  magnetic  field  Bq  and  the  mean  flow  field  Vq  are  assumed 
to  depend  on  the  coordinate  y  only  as  Bq  =  BQ^{y)e^  -\-  boe^^  and  Vq  =  Vo{y)ex,  where 
Boz{y)  =  By  ior  y  >  d,  Bo,{y)  =  Bq^  for  y  <  d,  Bo^{y)  for  \y\  <  d  \s  z.  given  arbitrary 
function,  6o  is  a  constant  and 

y/2 
Vyl2d 
I  -VI2. 
We  note  that  there  is  a  small  component  of  the  equilibrium  magnetic  field  6o  tangential  to 

the  flow  velocity  (6o  •C  .Boii)-  Here  this  system  has  been  transformed  in  such  a  way  that  the 
mean  velocity  Vo{y)  is  an  odd  function  and  the  flows  in  regions  (I)  and  (III)  are  given  by 
Vo/2  and  — Vo/2,  respectively.  Assuming  the  plasma  is  described  by  the  ideal  incompressible 
MHD  equations  (2.1a),  (2.1b),  (2. Id),  (2.1e),  and  V  •  v  =  0,  we  impose  a  small  perturbation 
whose  time,  x-  and  2-dependence  is  given  by  expi(u;i  —  k^x  —  k^z).  We  then  obtain  the 
following  set  of  equations: 

d_ 

'  dy 


ipijj'vij.  +  pviy^-vo    =    ikxPi  -iBQ2{k^Bii -kj.Bi^) ,  (3.19a) 


ipu'viy     =     —T-{pi  + Bq^Bu) -boi—T^ -\-ik:rBiyj 

-    ik,Bo,B^y,  (3.19b) 

ipu'vi^     =    ikzpi  +  Biy — r-^  +  ibo  {k^Bir  —  k^Bu) ,  (3.19c) 

dy 


dviy 


—  ikj:Vi,:  —  ik^Viz     =    0,  (3.19d) 


dBiy 


dy 

-  ik^Bi^  -  iKBu    =    0,  (3.19e) 


dy 
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iu'Biy    =     -i{boh  +  Bo,k,)viy,  (3.19f) 

iu.' Bu  +  viy—^     =     -i{bok^  +  Bo,k,)vu,  (3.19g) 

where  u>'  =  u  —  Vo{y)kj.  and  p  is  assumed  to  be  constant.  The  boundary  conditions  are 
that  all  perturbed  quantities  vanish  at  |y|  — >  oo.  We  now  consider  the  most  unstable  mode 
which  propagates  in  a  direction  perpendicular  to  the  field  (A:^  =  0)  in  the  region  (III)  or  the 
magnetosphere  [Southwood,  1968].  Using  k^  =  0,  we  obtain 

Biy    =    -^o,,,  (3.21) 

and 


_         ^    dB 

-Dlx       —       — 


ly 


kx     dy 

from  Eqs.  (3.19d),  (3.19f),  and  (3.19e),  respectively.   From  Eq.  (3.19a)  and  Eq.  (3.20),  we 

obtain 

.pu'  dviy       .   dvo 
kl     dy  dy 

Combining  Eqs.  (3.19b)  and  (3.21)-(3.23)  yields  the  following  equation  for  viy-. 


(pi  +  Bo5i,)  =  -i'^  -j^  -  ip-^v,y.  (3.23) 


)  [u''  -  kid)  +  2ktCl  (^) '  L:,  =  0,  (3.24) 


^  V  dy^ 

where  C^  =  bl/ p.  The  boundary  conditions  of  Eq.  (3.24)  are  given  by  i^iy  — >  0  as  |j/|  — v  oo. 
We  note  that  the  function  Boi(j/)  does  not  enter  Eq.  (3.24).  The  mean  magnetic  field 
component  Bqz  perpendicular  to  the  direction  of  the  wavenumber  vector  and  the  mean  flow 
velocity  does  not  affect  the  stability  in  the  incompressible  flow  case.    This  independence 
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of  the  stability  on  the  perpendicular  magnetic  field  is  also  known  for  the  case  of  a  zero- 
thickness  shear  layer  with  a  constant  perpendicular  magnetic  field  [Chandrasekhar,  1961]. 
Equation  (3.24)  can  be  significantly  simplified  with  the  use  of  new  variable 

Hy)  =  ^lyiy)/^'- 


With  this  new  variable  6,  Eq.  (3.19)  is  transformed  to 

d_ 

dy 

with  the  boundary  conditions  6{y)  — ♦  0  as  |y|  — >  oo.  To  simplify  Eq.  (3.25)  further,  we  also 
define  a  new  dependent  variable  ^  by 


I  {u;'\y)  -  kid)  1}  +  {u^'\y)  -  klCl)  6  =  0,  (3.25) 


(  =  kxy  -Q 

with  fi  =  Ljd/V.   In  the  regions  I  and  III,  where  u'  =  u  —  kxVo{y)  is  constant,  Eq.  (3.25) 
may  be  written  in  terms  of  ^  a^ 

^'<5  c  n 

The  solution  of  this  equation  in  the  region  I  (^  >  k^d  —  0.)  is  given  by 

6    =    Aiexp(-0         (^>M-n)  (3.26) 

and  the  solution  in  the  region  III  {(  <  —k^d  —  f2)  is  given  by 

8    =    >l3exp(0  {i<-k,d-^).  (3.27) 

Here  Ai  and  Az  are  constants.  In  the  region  II,  Eq.  (3.24)  becomes 


where 


a    = 


k    =     kxd 
A^    -     — 
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We  now  solve  Eq.  (3.28)  by  means  of  the  WKB  method.  For  the  range  of  the  variable  if  such 
that  1^^  —  a^  :7"^  0,  we  transform  (  to 


Za 


(-a 


^  +  Q 


and  therefore  we  have 


dt]    = 


d( 


-{e{r,)-a'y6  =  0. 


(3.29) 


Then  Eq.  (3.22)  is  transformed  to 

drj^ 

For  (f^  not  near  q^,  the  WKB  method  yields  the  asymptotic  solution  of  Eq.  (3.29)  in  terms 
of  ^  as 

6  =  A+  {e  -  a'y'^\xpi  +  A.  {e  -  a^y''' exp{-0.  (3.30) 

The  boundary  conditions  to  be  applied  at  z  =  ±(f  are  that  the  normal  velocity  viy{y)  and 
the  normal  component  of  the  stress  pi  +  BqBiz  be  continuous.  Since  u>'{y)  =  u  —  kxVo{y)  is 
continuous  a±  z  =  ±d,  the  continuous  normal  velocity  conditions  are 


[^].=±.  =  0, 


(3.31) 


where  [     ]^  indicates  the  change  in  the  value  of  the  enclosed  quantity  across  the  discontinuity 
surface  at  z.  Since  Eq.  (3.23)  may  be  written  as 

.  pu)'    dS 


(pi  +  BoBu)  =  -i- 


kx    dy' 


the  continuous  normal  stress  conditions  are  given  by 

=  0. 


d( 


(3.32) 

Equation  (3.32)  may  also  be  obtained  by  the  integration  of  Eq.  (3.28)  across  the  interfaces 
at  ,f  =  ±k^d-Q  =  ±k-9..  Equations  (3.31)  and  (3.32)  give  the  following  four  homogeneous 
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equations  for  Ai,  A+,  A_,  A^: 


Ai  +  A+e"  +  A-e-"     = 


Ai  +  A+ 


A3-A+ 


6 


1  1  e'^  +  ^_  I 
-We-'-A- 


0, 

A3-\-A+e-''  +  A.e''    =    0, 
+  l|e-'=     =    0, 


+  1     e'=     =    0. 


Here 


-4,    = 


'^/«^e-M„ 


(f)  -  )fc)y 


''^^^^e-^. 


-^3 

{n  +  k)v 

A 

=     e-"A, 

^ 

=     e-"^, 

6 

=  k-n, 

6 

=  -k-n. 

The  condition  that  the  equations  above  have  nontrivial  solution  Ai,  A+,  A-,  A3  gives  the 
dispersion  relation 


^"  +  02 


A:-2p(l  +  i^)-i(l-exp(-4fc)) 


+^'^^^"14^  -  ^'^  +  ^'  (^  -  ^'^P^-^^^))  =  0' 


(3.33) 


where  Q  =  2uid/V  and  A;  =  k^d.  Figure  7  shows  this  dispersion  relation  compared  to  the 
dispersion  relation  of  the  zero  thickness  magnetopause.  It  is  clear  that  the  growth  rate  of 
perturbations  with  large  k  is  significantly  reduced  by  the  effect  of  the  finite  thickness  of  the 
magnetopause  and  there  is  a  critical  kc  such  that  any  waves  with  k  >  kc  are  stable. 
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In  this  section,  we  have  reviewed  the  basic  theory  of  the  Kelvin- Helmholtz  instability 
at  the  magnetopause.  The  compressibihty  effect  on  this  instability  in  the  case  of  finite 
thickness  of  the  magnetopause  is  also  discussed  by  Ong  and  Roderick  [1972]  and  it  is  shown 
that  the  compressibihty  effect  does  not  alter  the  dispersion  relation  Eq.  (3.33)  significantly. 
The  required  condition,  Eq.  (3.17),  for  the  Kelvin-Helmholtz  instability  obtained  under  the 
assumption  of  zero  thickness  of  the  magnetopause  (which  gives  more  unstable  waves  than  the 
condition  obtained  under  the  assumption  of  finite  thickness),  however,  is  not  easy  to  satisfy 
[Hasegawa  and  Sato,  1989].  For  example,  with  Bi  =  2O7,  p2/(proton  mass)  =  0.1  cm~^  and 
(f>i  =  0,  then  Vi  >  2  X  lO^km/s,  which  is  rather  large  for  the  solar  wind  [Nishida,  1978]. 
Therefore  more  effort  is  needed  to  determine  a  significant  agent  of  wave  generation,  not 
only  considering  the  Kelvin-Helmholtz  instability  in  a  more  realistic  configuration  of  the 
magnetopause  but  also  taking  into  account  other  effects  such  as  dayside  reconnection. 

4.      Summary  and  Discussion 

In  this  review,  we  have  focused  on  Alfven  resonance  of  geomagnetic  field  lines  and  Kelvin- 
Helmholtz  instability  at  the  magnetopause  as  a  possible  mechanism  of  the  long-period  mag- 
netic pulsations  (pc3  ~  pc5).  In  Sec.  2,  we  discussed  the  theory  of  resonant  excitation  of 
the  Alfven  wave  at  a  particular  geomagnetic  field  line,  where  the  compressional  wave  from 
an  outside  source  is  converted  to  the  shear  Alfven  wave  due  to  the  nonuniformity  of  the 
magnetic  field.  In  Sec.  3,  based  on  the  incompressible  MHD  model,  the  dispersion  relations 
of  the  wave  on  the  magnetopause  are  derived  under  the  assumptions  of  the  zero  thickness 
and  the  finite  thickness  of  the  magnetopause.  It  is  shown  therein  that  finite  thickness  of 
the  layer  with  shear  of  the  magnetic  field  and  streaming  velocity  stabilizes  the  short  wave- 
length perturbations  which  are  unstable  in  the  zero  thickness  analysis.  In  these  theoretical 
treatments,  effects  of  magnetic  field  curvature  are  not  taken  into  account.  Despite  such 
ideahzation  of  the  models,  however,  it  is  beheved  that  essential  physical  mechanisms  of  the 
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magnetic  pulsations  are  contained  in  the  theories  presented  in  this  review. 

Apart  from  the  curvature  effects,  we  now  discuss  briefly  some  other  eff'ects,  particularly  on 
the  Kelvin-Helmholtz  instability.  Fejer  [1964]  and  Sen  [1964,  1965]  discussed  effects  of  com- 
pressibility on  the  Kelvin-Helmholtz  instability,  simplifying  the  problem  in  different  ways, 
and  drew  different  conclusions.  Fejer  finds  the  most  unstable  modes  propagate  parallel  to  the 
streaming  velocity  while  Sen  finds  the  most  unstable  modes  propagate  perpendicular  to  the 
unperturbed  magnetic  field.  The  relationship  of  these  conclusions  is  discussed  by  Southwood 
[1968]  in  a  more  general  treatment.  Both  Sen  and  Fejer,  however,  find  in  their  compressible 
MHD  models  that  there  is  an  upper  limit  of  the  streaming  velocity  for  each  normal  mode 
solution.  With  the  streaming  velocity  larger  than  this  upper  limit,  the  Kelvin-Helmholtz 
instability  is  stabilized.  This  upper  limit  of  the  streaming  velocity,  which  does  not  exist  in 
the  incompressible  model,  seems  to  have  httle  physical  significance  in  the  actual  stabihty 
problem  since  it  is  always  possible  to  choose  a  direction  of  the  wavenumber  in  such  a  way 
that  there  exists  an  unstable  mode  [Southwood,  1968].  Effects  of  more  realistic  configura- 
tions of  the  magnetosheath-magnetosphere  boundary  are  also  discussed  by  several  authors. 
For  example,  Lee,  Albano,  and  Kan  [1981]  studied  the  model  which  consists  of  three  uniform 
plasma  regions:  the  magnetosheath,  the  boundary  layer,  and  the  magnetosphere.  The  mid- 
dle layer  of  the  uniform  plasma,  as  a  matter  of  fact,  represents  a  boundary  layer  inside  and 
adjacent  to  the  magnetopause,  rather  than  the  magnetopause  itself.  This  boundary  layer 
consists  of  a  plasma  with  nearly  similar  temperature  and  flow  properties  to  the  solar  wind 
on  the  magnetosheath  due  to  permeabihty  of  the  magnetopause  [Haerendel  and  Paschmann, 
1982].  Lee,  Albano,  and  Kan  find  that  there  are  two  unstable  modes  in  the  magnetopause 
boundary  region:  one  is  excited  at  the  magnetosheath-boundary  layer  interface  (the  mag- 
netopause mode)  and  the  other  is  excited  at  the  inner  surface  of  the  boundary  layer  (the 
inner  mode).  While  the  magnetopause  mode  can  be  stabilized  by  the  magnetic  field  in  the 
magnetosheath,  the  inner  mode  is  found  to  be  unstable  most  of  the  time.  Figure  6  schemat- 
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ically  shows  the  inner  mode  at  the  interface  between  the  magnetosphere  and  the  boundary 
layer. 

Recently  Wolfe  et  al.  [1988]  discussed  the  penetration  of  the  magnetopause  mode  and 
the  inner  mode  deep  into  the  magnetosphere.  It  is  observed  that  the  monochromatic  surface 
waves  generated  at  the  magnetopause  have  sufficient  amphtude  deep  inside  the  magneto- 
sphere to  cause  a  field  fine  resonance.  As  first  pointed  out  by  Lanzerotti  et  al.  [1981],  the 
lack  of  understanding  of  the  mechanism  for  the  evanescent  waves  to  have  such  large  decay 
lengths  in  the  magnetosphere  was  a  major  difficulty  in  describing  the  long-period  magnetic 
pulsations  self- consistently.  It  is  found  by  Wolfe  et  al.  [1988],  however,  that  the  inner  mode 
may  have  its  energy  reduced  by  only  one  order  of  magnitude  over  a  penetration  depth  of 
about  SRe  under  some  conditions  while  the  magnetopause  mode  always  has  extremely  small 
decay  length  (~  IRe)  regardless  of  the  thickness  of  the  boundary  layer.  Thus  taking  into  ac- 
count the  finite  size  of  the  magnetosheath-magnetosphere  boundary  layer  and  the  associated 
inner  mode  may  yield  a  more  self-consistent  picture  of  the  long-period  magnetic  pulsations. 
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Name 

Frequency 

Name 

Period  or  Rise  Time 

SHF 

3—30  GHz 

Pc  1 

2«;...=0.2— S.Osec 

UHF 

0.3-3.0  GHz 

Pc2 

2n«<=5— lOsec 

VHF 

30-300  MHz 

Pc3 

2n'<>=IO— 45  sec 

HP 

3-30  MHz 

Pc4 

2«,<.j=45— I50sec 

MP 

0.3-3  0  MHz 

Pc5 

2jr/i'i=  150—600  sec 

LP 

30—300  kHz 

VLP 

3-30  kHz 

Pi  1 

T,=  1—40  sec 

ELP 

3-3000  Hz 

Pi2 

T,  =  40-150sec 

ULP 

<3  Hz 

sc.  si 

t,- 300  sec 

Classification  of  magnetospheric  signals  (after  Shultz  and  Lanzerotti,  1974). 
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Classification  of  magnetic  pulsations  (after  Saito  et  al.,  1976). 


TABLE   2 


2« 


ALTITUDE  ■ 


■km  L  -    V  AL  U  E 


soo  Doo  noo   4000    2         >     *  ?        '9 


Pes 


Pe4 


PcS 


Pe2 


"10'       »o^^       »o^       •?       »o^       *^       *** 
6E0MA0NETIC    LATITUDE 


Fundamental  period  T  of  the  magnetic  field  oscillation  given  by  Eq.  (1.1).  The  dipole 
magnetic  field  is  assumed,  where  field  lines  are  labeled  by  the  geomagnetic  latitude  of 
their  intersection  with  the  earth  or  by  the  equatorial  crossing  distance  (L- value)  (after 
Saito  et  al.,  1976). 

FIGURE   1 
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m<0 


m>0 


Tilt  of  the  major  axis  near  the  resonant  field  Une  at  the  equatorial  plane  (after 
Hasegawa  and  Chen,  1974). 


FIGURE  2 
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Schematic  diagram  of  the  wave  amplitude  and  sense  of  polarization  vs.  radius  in  the 
equatorial  plane  for  the  local  morning.  When  projected  along  the  field  lines  to  the 
ground,  it  can  be  viewed  as  a  function  of  the  latitude.  C  and  E  stand  for  circular 
and  elliptical  polarization,  respectively,  and  Uq  is  the  location  of  the  resonant  field  line 
(after  Hasegawa  and  Chen,  1974). 

FIGURE  3 
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The  diurnal  variation  of  the  sense  of  polarization  for  quasi-monochromatic  pulsation 
(frequency  ~  5mHz)  (after  Samson  et  al.,  1971). 


FIGURE  4 
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aV 


(I)  Magnetosheoth 


(I)  Magnetosphere 


■*►  X 


(a) 


(b) 


Coordinate  system  of  the  problem  (a)  the  unperturbed  magnetopause  is  the  x-z  plane, 
(b)  the  definition  of  the  angles  <^  and  V  (modified  after  Sen,  1963). 


FIGURE   5 
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Shear  layer  model  (modified  after  Ong  and  Roderick,  1972). 


FIGURE  6 
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A— 00  (b.>0) 


Growth  rate  vs.  wavenumber  in  the  case  where  finite  thickness  of  the  magnetopause  is 
taken  into  account  (after  Ong  and  Roderick,  1972). 


FIGURE   7 
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Schematic  interpretation  of  the  inner  mode.    The  shaded  area  is  the  boundary  layer 
and  its  wavy  inner  edge  is  indicated  by  the  dotted  Une  (after  Hones  et  al.,  1981). 


FIGURE   8 
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